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. Abstract: We review how to determine \Vub\ from inclusive semileptonic B decay using 1 < ^ 1 

, combined cuts on the Icptonic and hadronic invariant masses to ehminate the 6 — > c 

I background. This leads to a determination of \Vub\ with theoretical uncertainty at the 

5-10% level. 
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The magnitude of Vub determines the length of one of the side of the unitarity triangle I — ^ 
' and is therefore of great importance to overconstrain this triangle. The inclusive decay 

rate B Xulv is directly proportional to \Vub'\^ and can be calculated reliably and with 
^ . small uncertainties using the operators product expansion (OPE). Unfortunately, the ~100 

times background from B — > Xclv makes the measurement of the totally inclusive an almost 
impossible task^. Several cuts have been proposed in order to reject the 6 — >■ c background, 
however care has to be taken to ensure that the decay rate in the restricted region of phase 
space can still be predicted reliably theoretically. The cut which is easiest to implement 
experimentally is on the energy of the charged lepton, requiring E^ > {rn?^ — rnj^) /2mB- 
Unfortunately, this cut restricts the remaining region of phase space too much for the OPE 
to still be valid. Instead, a twist expansion has to be performed [||, and at leading order 
the decay rate is determined by the light cone distribution function of the B meson, with 
subleading twist corrections suppressed by powers of l/mb [Q]. This distribution function 
can not be calculated perturbatively and has to be determined experimentally, for example 
from the photon energy spectrum in — > Xg-^ [^. The same light cone distribution 
function also determines the rate in the presence of a cut on the hadronic invariant mass 
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Q, which has the advantage of keeping most of the b ^ u events. In this talk I will show 
that using a cut on the leptonic invariant mass results in a regions of phase space that is 
free from b ^ c background and which can be calculated reliably in the standard OPE. 

Using a cut on the leptonic invariant mass = {pe + Pu)^ to measure iKibP was first 
proposed in Q, and it was shown that requiring > — leaves a region of phase 
space free from b ^ c events, which can be calculated using the usual OPE. The number 
of events surviving such a cut on were calculated in ^ and depending on the exact 
value of the cut chosen, the fraction of events surviving the cut is between 10 and 20%, 
with uncertainties on \Vub\ ranging from 15% for q^^^^ = m\ — mj^ = ll.GGeV^ to 25% for 

= 14GeV2. 

From Fig. |^ one can see that the q"^ cut keeps a 
subset of the phase space available to a cut on the 
hadronic invariant mass mx- Thus, a cut on mx 
keeps a much larger fraction of the events, up to 
80%. However, the fact that the former cut avoids 
the region of low mx and low q^ where the structure 
function dominates, is why the decay rate in the 
presence of a cut can be calculated while in the 
presence of a pure mx cut it can not. However, it 
was observed in |^ that combining the cuts on q^ 
and on mx allows for a much larger fraction of events compared to a pure cut, while 
keeping the partial rate calculable using the local OPE. 

The integrated rate with a lower cut g^^^ on q'^ and an upper cut mcut on mx may be 
written as 

r _ G^|Kfe|^(4.7GeV)5 

/o ^ dq^ds ^ 1927r3 



(GeVl 




Figure 1: The dalitz plot for and sh- 



df 



(1) 



where where q = q/m},, s = [v — g)^ is the rescaled partonic invariant mass, v is the 
four-velocity of the decaying B meson, and 
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for mcut > ms - mb^/q^ , 
for ml^,<{mB-m,e){mB-m,), (2) 

1 I otherwise . 



The hadronic invariant mass mx is related to q^ and s by 



m\ 



smBmi, + {ms - mh){mB 



(3) 



G{qcui,mcut) is the ratio of the semileptonic b —> u width with cuts on and mx to the 
full width at tree level with mb = 4.7 GeV. The fraction of semileptonic b ^ u events 
included in the cut rate is ~ 1.21 G(gcut) "^cut)- Note that the m^ prefactor, a large source 
of uncertainty, is included in G{q'^^^,mcut)- The theoretical uncertainty in \ Vub\ is therefore 
half the uncertainty in the prediction for G(gcuti "^cut)• 
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In Ref. @ this function G{q^^^,mcut) was cal- 
culated up to corrections of order terms not 
enhanced by /3o, order [^qcti / fnhf' terms propor- 
tional to derivatives of 5{s), and higher order terms 
in both series. The result for the function , mcut 

is shown in Fig. ^ for three different values of mcut 
as a function of g^ut- 

The uncertainties in the OPE prediction for 
G((7cut' ""^cut) come from three separate sources: per- 
turbative uncertainties from the unknown full two- 
loop result, uncertainties in the h quark mass and 
uncertainties due to unknown matrix elements of 
local operators at 0(l/m^) in the OPE. In the fol- 
lowing I will consider each of these uncertainties separately as the fractional errors on 
G{q^^^,mcut)- The fractional uncertainty in \Vuh\ then is one half of the resulting value. 
Starting with the perturbative uncertainties, we may estimate the error in the perturba- 
tion series in two ways: (a) as the same size as the last term computed, the order a^/3o 
term, or (b) as the change in the perturbation series by varying ^ over some reasonable 
range. These are illustrated in Fig. |^ (a) and (b), respectively. In Fig. |3|(b) we vary the 
renormalization scale between /i = 4.7 GeV and ^ = mfe/3 ~ 1.6 GeV, and plot the change 
in the perturbative result. For a given set of (/^ut aiid nicut, we take the perturbative error 
to be the larger of (a) and (b). 



Figure 2: G((7^„t, mcut) as a function 
of g^uti for TOcut = 1-86 GeV (solid line), 
1.7 GeV (short dashed line) and 1.5 GeV 
(long dashed line). 
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Figure 3: (a) The 0(e) and ©(e^lm) contributions to G{q^^^.,m cut) (normalized to the tree level 
result) for Wcut = 1.86 GeV (solid lines), 1.7 GeV (short dashed lines) and 1.5 GeV (long dashed 
lines), (b) Scale variation of the perturbative corrections: The difference between the perturbative 
corrections to G{qcut, nicut), normalized to the tree level result, for /i = 4.7 GeV and /i = 1.6 GeV. 



The partially integrated rate depends sensitively on the value of the b quark mass due 
both to the factor in G{q'^^^,mcut) and the cut on g^, as stressed in Currently, the 
smallest error of the 15" mass is quoted from sum rules ^, 

value m^^ = 4.69 it 0.03 GeV by fitting an optimized linear combination of moments of the 
e'^e~ —>■ bb spectrum, which may underestimate the theoretical error |^]; the authors of 
Q cite a similar central value with a more conservative error of ±0.08 GeV. In Fig. | we 
show the effects of a itSOMeV and a itSOMeV uncertainty in on G{q'^^^^,mcut), using 



10 1 . Ref. obtains the 
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the central value = 4.7 GeV. 
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Figure 4: The fractional effect of a ±80MeV and ±30MeV uncertainty in ml^ on G{q^.^^, nicut) 
for rricut = 1.86 GeV (solid line), 1.7 GeV (short dashed line) and 1.5 GeV (long dashed line). 

Finally, there are uncertainties from unknown l/mf contributions to mcut)- We 

distinguish two contributions, the first originating from the local OPE and the second 
from weak annihilation (WA) [|l2|. The latter terms vanish in the factorization limit, 
however they are enhanced by a factor of IGtt^. The two contributions are shown in Fig. |5|, 
assuming a 10% violation of factorization, with the uncertainty from WA dominating. Since 
the WA contribute in general differently to charged and neutral B decays, an indication of 
the size of the WA effect can be obtained by comparing the measurement of mcut) 
from charged and neutral S's. 
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Figure 5: Estimate of the uncertainties due 
to dimension-six terms in the OPE as a func- 
tion of q^^^ from weak annihilation (WA) 
(solid line) and other operators (dashed line). 
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Figure 6: Effect of a model structure func- 
tion on G'((j^„t, mcut) as a function of for 
m^^t = 1.86 GeV (solid line), 1.7 GeV (short 
dashed line) and 1.5 GeV (long dashed line). 
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Since the structure function becomes more important as the cut on is lowered, we 
still need to estimate the effect of the distribution function to determine how low may 
be decreased. To leading twist, this is obtained by smearing the h quark decay rate with 
the distribution function f{k^). The best way to determine f{k+) is from the B Xg^ 
photon spectrum, however in the absence of precise data, we will use a simple model to 
estimate the effects of the structure function. In Fig. ^ we plot in such a model the effect 
of the structure function on G(g^^^,mcut) as a function of q^^^, for three different values of 

"Icut- 

Using these results the strategy to extract \Vub\ in a model independent way is: 
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• make the cut on mx as large as possible, keeping the background from b ^ c small 

• for a given cut on mx, reduce the cut as low as possible, keeping the contribution 
from the b quark structure function, as well as the perturbative uncertainties, small 

In Table ^ we show the results 
for two different combinations of 
cuts. From this table one can see 



(9?ut>"lcut) 

(GcV2, GcV) 


Fraction 
of Events 


A/Kb 


AKb 
±80/30 


Combined cuts 


(6, 1.86) 
(8,1.7) 


46% 
33% 


-2% 
-3% 


8%/5% 
9%/6% 


Pure cut 


(11.6, 1.86) 


17% 




15%/12% 



Table 1: Fraction of events, effect of structure function 
and uncertainty on \Vub\ for two different combination of 
cuts, in comparison with a pure cut. 



that the number of events included 
using the combined cuts is about 
a factor of three larger than for a 
pure cut and that the uncertain- 
ties on \Vub\ can be reduced by a 
factor of two. 

In summary, combining a cut 
on with a cut on mx allows for 

a model independent determination of \Vub\ using up to 45% of the total b ^ u events 
with uncertainties at the 5-10% level. This is a significant improvement over other meth- 
ods presented to date, which either have strong model dependence or larger theoretical 
uncertainties. 
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